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Abstract. We introduce a transport approach at fixed shear viscosity to entropy ratio rj/s 
to study the generation of collective flows in ultra-relativistic heavy-ion collisions. Transport 
theory supplies a covariant approach valid also at large rj/s and at intermediate transverse 
momentum pr, where deviations from equilibrium is no longer negligible. Such an approach 
shows that at RHIC energies a temperature dependent rj/s enhances significantly the V4/V2 
respect to the case of constant rj/s. Furthermore if NJL chiral dynamics is self-consistently 
implemented we show that it does not modify the relation between V2 and rj/s. 



1. Introduction 

The RHIC program at BNL has shown that the azimuthal asymmetry in momentum space, 
namely the elliptic flow V2, is the largest ever seen in HIC suggesting that an almost perfect 
fluid with a very small shear viscosity to entropy density ratio, rj/s, has been created [U [2]. 
Nonetheless it has been found that dissipative effects cannot be neglected and even a small shear 
viscosity to entropy ratio rj/s produce sizeable effect increasing with the transverse momentum 
Pt of the particles [3]. This has triggered a lot of activity in developing a relativistic theory 
of viscous hydrodynamics. Viscous corrections to ideal hydrodynamics are indeed large and a 
simple relativistic extension of first order Navier Stokes equations is affected by causality and 
stability pathologies [5]. It is therefore necessary to go to second order gradient expansion, and 
in particular the Israel-Stewart theory has been implemented to simulate the RHIC collisions 
providing an upper bound for r//s < 0.4 [JJ. Such an approach, apart from the limitation to 
2+1D simulations, has the more fundamental problem of a limited range of validity in rj/s and 
in the transverse momentum px- 

On the contrary a relativistic transport approach has the advantage to be a 3+ ID approach 
not based on a gradient expansion in viscosity that is valid also for large r//s and for out of 
equilibrium momentum distribution allowing a reliable description also of the intermediate pt 
range where the important property of quark number scaling (QNS) of V2{pt) has been observed 
|10j . In this pt region viscous hydrodynamics breaks its validity because the relative deviation 
of the equilibrium distribution function 5f/f eg increases with p T becoming large already at 
Pt > 3T ~ lGeV [4J. Furthermore Boltzmann-Vlasov transport theory distinguishes between 
the short range interaction associated to collisions and long range interaction associated to the 
field interaction, responsible for the change of the Equation of State (EoS) respect to that of a 
free gas. This last feature allows to unify two main ingredients that are relevant for the formation 
of collective flow. In ideal hydrodynamics the V2(pt) depends strongly on the EoS namely on 



the sound velocity (? s = dP/de [8], while the mean free path A is assumed to be vanishing. In 
the parton cascade approach the EoS is fixed to be the one of a free gas c 2 s = 1/3 = P/e, on the 
other hand the mean free path A = 1 /pa is finite. In the first stage of RHIC the two different 
approaches were able to account for the large t>2(pr) observed; in particular the parton cascade 
with large scattering cross section predicted the saturation of V2(pr) vs px [9]. Anyway once 
viscosity is finite both a finite A and the EoS are important for the generation of the momentum 
anisotropies. 

The basic equation of transport for the (anti-) quark phase-space distribution function 
under a scalar field interaction can be written as: 

p^/^s.p) + M(x)d^M(x)d^f ± (x,p) = C(x,p) (1) 

where the first term is related to the free streaming, the second term represents the effect of a 
scalar field modifying the e = 3P relation (giving a finite interaction measure, see Fig. 2 (left)) 
and the last term is the effect of the collisions directly associated to a finite A and therefore to 
a finite rj/s. 

In the following we will show two results obtained with such an approach: one related with 
the large v±/v\ ratio observed at RHIC and the other discussing the role of the NJL chiral 
dynamics on the relation between vi(jpr) an d rj/s. We will show calculations for Au + Au at 
200 AGeV; the density distribution in coordinate space is given by the standard Glauber model. 
The maximum initial temperature is T = 340 MeV and the initial time is To = 0.6 fm/c as 
usually done also in hydro dynamical calculations. 

2. The large v^/v\ ratio 

We first consider the reduction of the transport approach to a cascade, i.e. neglect the field 
interaction, included in the next section. In such a case the mass M(x) = in Eq.(l). It has 
been shown that a fluid at finite rj/s can be studied within a transport approach renormalizing 
locally the scattering cross section according to r\j s-ot r =< V > /15p, for details see Ref.s [BJEIE]. 
The system of course expands and cools and when the energy density is below e = IGeVj fm 3 
it starts the cross-over region and then the hadronic stage. During this phase we increase the 
rj/s of the medium according to Fig. 1 (left) in order to have finally in the hadronic phase an 
rj/s = 8/4tt ~ 0.65 in agreement with several calculation of rj/s. 

At RHIC it has been observed v^ipx) /v^ipr) ratio of about 0.8-0.9 quite larger than first 
ideal hydro dynamical predictions of about 0.5. More realistic hydro dynamical calculations show 
that the ratio can be even lower than 0.5 and viscosity does not strongly affect such a ratio but 
can bring it up to about 0.6. Furthermore in Ref.[TT] it has been shown that such a ratio is 
almost independent on the initial eccentricity, on the viscosity, on the particle species. Therefore 
even if v± and v-i can depend on several details such a ratio is a quite stable one and its large 
value should be considered as an oppurtunity to have some new information on the system not 
yet spot. In Ref.[7| we have shown that indeed the formation time of is larger than the v% 
one and therefore it is more affected by dynamics of the later stage of the reaction. Thanks 
to the difference in the formation time we argue that such a ratio is significantly affected by 
the temperature dependence of the rj/s. Therefore we have performed two calculations one 
with a constant Airrj/ s = 1 in the QGP phase (dashed line) and the other, shown by solid line 
in Fig. 1 (left), with a T or e energy density dependence inrj/s = T/Tq = (e/eo) 1//4 where 
eo = 1.7 GeV/ fm 3 gives the beginning of the cross over region. We remind that an increase of 
rj/s is expected moving away from the phase transition region on general considerations |12j . 

In Fig.l (left) we show the results of the two calculations for Au + Au at 200 AGeV. For a 
constant rj/s we find a value of about 0.6 while in the case of a temperature dependent rj/s, it 
increases getting very close to experimental data, v±/v\ ~ 0.8. The reason for the sensitivity of 
this ratio on the T dependence is that the elliptic flow t>2 mainly develops earlier respect to V4 
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Figure 1. Left: Dependence of rj/s on the energy density e implemented in the transport code. 
Right: v^/v\ for Au + Au collisions in the mid-rapidity region \y\ < 1 at b= 8 fm; dashed line 
is the result for a constant Airrj/s = 1, the solid line is the result for a temperature dependence 
of the rj/s like shown by solid line on the left figure. 



and therefore on the average at larger rj/s (for the solid line case), while V4 develops more when 
the 77/s is smaller hence is less suppressed by dissipative effects. 

3. Impact of NJL Dynamics 

To start the investigation of the role of the field interaction we have implemented the NJL 
dynamics that introduces the generation of a finite mass associated to chiral symmetry breaking. 
Therefore particles do not have a vanishing mass but a finite one according to the mass gap 
equation of the NJL model: 

M(x)-m f d 3 p 1 -f~(x,p) - f + (x,p) 

that determines the local mass M{x) at the space-time point x in terms of the distribution 
functions f^(x,p). For a gas at finite mass e — 3P is no longer vanishing as shown in Fig. 2 
(right) for several NJL parametrization. 

Eqs.(Q]) and ^) form a closed system of equations constituting the Boltzmann-Vlasov 
equation associated to the NJL Lagrangian. The space-time dependence of the mass M a {r,t) = 
m — 2g(ipip) influences the momenta of the particles because the finite gradient of the condensate 
generates a force which changes the momentum of a particle proportionally to V r {ipip). The 
last is negative because the phase transition occurs earlier in the surface of the expanding QGP 
fireball. Therefore the phase transition which take place locally results in a negative contribution 
to the particle momenta that makes the system more sticky respect to a free massless gas. To 
study a fluid at finite rj/s and NJL field interaction one has to extend the simple relation 
ftr " v/ s = < P > /15p valid for a massless gas to the general case of massive relativistic particles 
and this has been discussed in Ref . [13] . 

In Figf21 we show the elliptic flow as a function of pr- We find that the presence of an 
NJL-field that drives the chiral phase transition suppress the V2{pt) by about 20% at px > 1 
GeV. This would imply the need of a parton scattering cross section atr even larger than that 
estimated with the cascade model which was already quite larger than the pQCD estimates [S] 
in order to describe the data. On the other hand the mean field modifies both the local entropy 
density, reduced by the mass generation, and the shear viscosity that increases respect to the 
massless case. Considering that one of the main goal is to determine the rj/s of the QGP, we 



have investigated what is the action of the mean field once the rj/s of the system is fixed to be 
the same with and without the NJL field. This of course implies different values of the cross 
section in the two cases. 




Figure 2. Left: Interaction measure shown as a function of temperature T for three different 
NJL parameter sets. Right: Average elliptic flow as a function of time for Au + Au collisions in 
the mid-rapidity region \y\ < 1 at b= 7 fm; see text. 

The results for 4irr]/s = 1 are shown by dashed lines in Fig. [2] (right), the (black) dashed 
line is the case with only the collision term (cascade) while the (green) dashed line is the case 
with the field. We can see that once the rj/s is fixed there is essentially no difference in the 
calculations with and without a field dynamics included. This is a key result that shows that 
even in a microscopic approach that distinguishes between the mean field and the collisional 
dynamics the V2(pt) is mainly driven by the rj/s of the fluid. In other words we have found 
that in a microscopic approach the rj/s is the pertinent parameter and the language of viscous 
hydrodynamics is appropriate. Of course this does not mean that V2(pt) in the transport theory 
is the same as in viscous hydrodynamics (especially at px > 1-5 GeV), but that the direct 
relation between V2(pt) an d f]/s is not modified by the NJL field dynamics. This confirms the 
validity of the studies pursued till now even if they miss an explicit mean field dynamics and/or 
the chiral phase transition. 
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